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Abstract 

We analyze the structure of the periodic trajectories of the K-system gen- 
erator of pseudorandom numbers on a rational sublattice which coincides with 
the Galois field GF[p\. The period of the trajectories increases as a function 
of the lattice size p and the dimension of the K-matrix d. We emphasize the 
connection of this approach with the one which is based on primitive matrices 
over Galois fields. 
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1 Introduction 



Nowadays the Monte-Carlo method has a wide range of applications and the quality 
of pseudorandom numbers being used plays an important role. Different principles 
and algorithms have been suggested in the literature to generate pseoudorandom 
numbers and to check their properties Q. The development of the ergodic theory 
[0. |3|, ^, [J and the progress in understanding of nonlinear phenomena together with 
the increasing power of modern computers open a new era for applications |7], |8|, |9j . 

In the articles ]TD| the authors suggested to use many-dimensional Kolmogorov 
K-systems to generate pseudorandom numbers of high quality. K-systems are the 
most stochastic dynamical systems, with nonzero Kolmogorov entropy and their 
trajectories are exponentially unstable and uniformly fill the phase space 0, |3], ||, [L^, 
|TB|, PH- It was suggested to use the coordinates of these trajectories as a sequence of 



pseudorandom numbers [10]. From this point of view the most successful inversive 



congruential generator Jy] used so far can be considered as a one-dimensional K- 



system and it was pointed out that this fact explains its exceptional properties flO" 

For the application of this idea it is important to have such K-systems for which 
the phase space is limited by a unit d- dimensional torus, because in that case the 
coordinates of the trajectories can be used directly without any additional trans- 
formations. Two types of K-systems have been suggested for these purposes: toral 



automorphisms [[LQ, 11] and many-dimensional Sinai billiard which is defined inside 



a unit d- dimensional torus |TJ]. 

In the case of toral automorphisms a unit d- dimensional torus U d plays the role of 
a phase space and the K-system is represented by a rf-dimensional matrix - K-matrix 
- which acts on the vectors from U d generating trajectories uniformly distributed 
over the torus Il d . The coordinates of these trajectories are used for Monte-Caro 
simulations [|1(], |TTJ. The properties of this new class of matrix generators were 
investigated by different criterion including Kolmogorov discrepancy . In all 
cases it shows good statistical properties [TT[ . 



The aim of this article is to estimate the period of the trajectories which are 
used to produce pseudorandom numbers generated by a K-system. It is clear, that 
only periodic trajectories of K-systems can be simulated on a computer, because 
trajectories on a computer are always on a finite rational sublattice Z d of the phase 
space U d . Thus we have to consider the system on rational sublattice Z d of a unit d- 
dimensional torus and particularly on sublattices with prime basis p p3| , [2~4| , |25 , 27 , 
EE |36| , |32| , pHJ . These sublattices are equivalent to Galois fields GF[p] and all four 



elementary arithmetical operations can be carried out unrestrictevely pi], |20, |22 



Analyzing trajectories of a K-system on a Galois sublattice GF\p] one can see 
that in order to have trajectories with large period K-matrix should have an eigen- 
values in high extensions GF[^/p] of the field (notation used in mathematical litera- 
ture is GF[p d }). This property makes them very close to so called primitive matrices 
which have been considered by Niedereiter [16|, |Tj]] to generate pseoudorandom num- 
bers. We refer to the book of Niederreiter |19[ and to the survey article |18| for recent 
references. The main idea of his approach is to use a primitive matrices on a given 
Galois field GF[j¥p] to generate pseudorandom numbers of very large periods. This 
approach guarantees the large period of the series. In addition the fascinating result 
of Niedereiter [[13] allows to estimate the uniformity of maximally long trajectories 



in terms of Kolmogorov discrepancy D^. 

Thus these two approaches are very close to each other on Galois sublattices 
and the main question is: whether one can have the matrices with both properties 
at the same time? The determinant of a K-matrix should be equal to one while 
the determinant of a primitive matrix is different from one, thus these properties are 
incompatible. The main point, which we would like to stress here, is that nevertheless 
one can construct K-matrices which have a primitive matrix as submatrices. In that 
case the trajectories are still very long as in the case of primitive matrices, but at 
the expense of appearance of trajectories with short period. Excluding them from 
initial data we guarantee that the trajectories are maximally long and at the same 
time belong to a K-system. We suggest specific matrices with these properties which 
can be used for practical simulations. 



2 Trajectories of K-system on a rational sublat- 
tice 

Let us pass to the details of the algorithm. The matrix generator is defined as 
[0,0, 

X (n+i) = A . X ( n \ (modi), (1) 

where A is d x d dimensional matrix with integer matrix elements a^j and determi- 
nant equal to one 

Det A = l, (2) 

and X® = (X{ 0) ,...,xi 0) ) is an initial real vector. The last condition provides 
phase-space volume conservation. The automorphism (1) forms the .fT-system of 
Anosov if and only if all eigenvalues of the matrix A are in modulus different from 
unity |, 1, | 

I Ai| ^1, i = l, -,d (3) 
The trajectory of the K-system (1) 



Xo, Xi, X 2 .... 

represents the desired sequence of the pseudorandom numbers flit)] . 

This approach allows a large freedom in choosing the matrices A for the K-system 
generators and the initial vectors QlOj . Specific choices suggested in [ffOl O, are 



/ 2,3,4, ,d ,1 \ 

1,2,3, ,d- 1,1 

1,1,2 d-2,1 



1,1,1, ...,2,3, 4,1 
1,1,1, ...,1,2,2,1 
1,1,1, ...,1,1,2,1 
1,1,1, ..,1,1, 1,1 ) 



A n 



( 0, 1 , 
0, , 1 



0, 



\ (-1) 



, o ,.. 
, ai, ct2, 








, 1 



(4) 



The first matrix has the advantage to be well defined in any dimension d and it has 



a very large Kolmogorov entropy [10| which is given by the Anosov-Sinai formula 



h(A d 



ln ^k- 

|A te [>1 



The entropy h defines the number 7r(r) of the periodic trajectories with period less 



or equal to r |30|, |31], p2 



tt(t) 



e hr 



when r — > oo, thus the number of "available" trajectories increases with entropy. 
The second one has a very simple expression for its characteristic polynomial 



A — ad-i A 



d-l 



ai A + (-r 



and for its eigenvalues Ai,...,A<2 we have Ai • ■ • A^ = 1,.., Ai + ... + A^ = a^-i- 
These formulas allow to choose eigenvalues and then to construct K-matrices. This 
correspondence between matrices and polynomials has a wide range of applications 
in algebra and number theory [p0[| . in the given case DetAd = 1 to fulfil K-condition 

(2) m. 

Let us consider trajectories of the system (1) with an initial vector which 
has rational coordinates |23|, E4I |25L E71 EM Wl 



Pi P2 Pd 



Qd 



(5) 



It is easy to see, that all these trajectories are periodic orbits of the Anosov map 
(1), because matrix elements integer. Indeed, if we consider the sublattice 

of unit torus U d with rational coordinates of the form q/p where p is the least 
common multiple of p's 



X 



Si 
p 



Si 
p 



p 



0<qi<p-l 



then the multiplication, summation and (mod) operations (1) will leave the trajectory 
on the same sublattice. The total number of vertices on this sublattice Z d is 



(total number of verteces) = p d , 



therefore the period t p of the trajectories on = Z p <g> ... ® Z p , where Z p = 
{0, 1, ...,p — 1} is always less than p d 



Thus the periodic trajectories of this system (1) with the initial vector (5) coincide 
with a subset of the points of rational sublattice Z^ and our goal is to find conditions 
under which the period of the K-system will be as large as possible. 

Let us show that on every given sublattice Z^ Anosov map (1) reduces to (mod 
p) arithmetic. Indeed on sublattice Z^ the Anosov map A (1) can be written as 

q (n+l) (n) 

= y^ciu - J — , {mod 1) 

p j P 

and is equivalent to (mod p) arithmetic on the lattice with integer coordinates q^ 
which are in the interval [0,p — 1] 

q(n+i) = J2a id ql n \ (mod p). 

j 

Thus the images of the periodic trajectories on a unit torus U d appear as trajectories 
on the integer sublattice Z d and all operations can be understood (mod p). The 
most important thing is that now all operations become commutative. 

To estimate the period of the trajectories on a rational sublattice it is essential 
to consider those sublattices for which p is the prime number, we mean that p\ = 
... = p d = p [E3, E4l ER E7L E6|, p2l E8|. In that case the integer sublattice gains an 



additional structure and becomes the Galois field GF[p] and all operations reduce 
to arithmetic ones on Galois field. The benefit to work on Galois field is that four 



arithmetic operations are well defined on that sublattice |H[ . 

In this way we can consider every coordinate qi , i = l,...,d as belonging to 
Galois field GF[p] = {0, 1, p — 1}, where p is a prime number and consider the 
sublattice direct product of Galois fields. 

Z$ = GF\p]®...®GF\p]. 

As we already mentioned, this reduction of a dynamical system (1) to a dynam- 
ical system for which the Galois field plays the role of the phase space makes all 
operations commutative in the sense that 

{A{A X}} = {A 2 X}, 

where {...} means mod operation. The commutativity of the multiplication and 
(mod) operation on the Galois sublattice means that the periodic trajectory 

{•M-i {A .v }...}} ,v 

can be represented in the form 



{A Tp X} = X. 



This equation allows to understand the relation between eigenvalues of the matrix 
A and the period of the trajectories. Indeed let us consider the eigenvalue problem 
for the matrix A on a Galois sublattice 

A X = X X, (6) 

then the period of the given trajectory r p can be understood as a degree of power 
on which the A reduces to identity (mod p) 

\ Tp = 1 (mod p). (7) 

The period of the trajectory on a Galois sublattice GF[p] is equal therefore to the 
power t p in which the eigenvalue of the matrix A reduces to identity. It is obvious 
that the same matrix A will have different periods on different Galois fields and that 
this period depends on the given prime number p,the dimension of matrices d and 
the initial vector X . 

3 Eigenvalues of the generator and the period of 
the trajectories 

Thus the actual value of the period r p naturally depends on the form of eigenvalues 
A and of the prime number p. Here we can distinguish different cases: 

i) . The eigenvalue A coincides with one of the elements of the Galois field GF[p\. 
In that case the period r p depends on the fact whether eigenvalue coincides with a 
primitive element of the Galois field. All elements of the field GF\p\ can be con- 
structed as powers of primitive element g and g v ~ x = 1. If the eigenvalue coincides 
with the primitive element of the Galois field , 

A = g, where g is a primitive element of GF[p], (8) 
then the period is maximal and is equal to r p = p — 1 

A p_1 = 1, (mod p). (9) 

Therefore to get the maximal period in the case i) one should have an eigenvalue 
equal to the primitive element g. If A does not coincide with the primitive element 
g, then the period is simply smaller and is equal to (p — l)/m where m is a divisor 
of p — 1 . 

ii) . The eigenvalue does not coincide with any element of the Galois field GF[p]. 
This may happen because Galois field is arithmetically complete, but it is not al- 
gebraically complete, therefore one can have the situation when the solution of the 
characteristic polynomial of the K-matrix is not in the field GF[p\. In that case one 
should ask, whether it is an element of the quadratic extension GF[^p\ or of higher 
extensions. The quadratic extension of the Galois field consists of the numbers of 
the form a+b^Tg where a, b are the elements of field GF[p], g is the primitive element 



of GF\p] and ^Jg is a square- free integer. The primitive element of the GF[<yp] has 
the period equal to p 2 — 1 ]2l| . 



Thus if the eigenvalue is an element of the quadratic extension and coincides 
with its primitive element h 



A = h, where h — hi + h 2 - s fg is a primitive element of GF[y/p\, (10) 
then the period is equal to r p = p 2 — 1 

X p2 ^ = 1, {mod p). (11) 

iii). In general the characteristic polynomial of the K-matrix is of order d and the 
eigenvalue may belong to high extensions GF[^/p] of the Galois field. The elements 
of GF[j¥p\ have the form a + bh + ... + eh 1 where a, b, e are the elements of 
GF[p] and h is a primitive element of GF[^/p] ]20|, [H], p2] . If the eigenvalue A 
coincides with this primitive element 



A = h where h is a primitive element of GF[tfp\, (12) 



then the period is equal to r p = p d — 1 EO, |2Tj, |22 



\ p - 1 = 1, (mod p). (13) 

This analysis demonstrates an important fact that in order to have a large period 
on a sublattice GF[p] one should have K-matrices with eigenvalues in high extensions 
of the field. 

4 Generators with largest period 

In the previous sections we described the trajectories of the K-system on the rational 
sublattice Z d and particularly on a Galois field, that is when p is a prime number. 
We have seen that the period of the trajectories depends on the "order" of the 
corresponding eigenvalue and the period is as large as the order of the extension 
of the field to which belongs the eigenvalue. The question is: can we construct a 
K-matrices with eigenvalues in high extensions of the Galois field and how many of 
them can simultaneously belong to a maximal extension GF[tfp] ? 

We should remark that the d-dimensional matrices A with all eigenvalues in 
GF[^/p] are well known in number theory and correspond to so called primitive 
matrices of the field GF[^/p] and the determinant of primitive matrices is not equal 



to one |20| , |21| , |22| . Therefore the K-matrices which have the determinant equal to 
one can not coincide with the primitive matrices, but as we will see one can construct 
K-matrices with d — 1 eigenvalues in GF[tfp] and only one in GF[p\. This means 
that most of the trajectories will have the maximal period r p = p d — 1 and only few 
of them (corresponding to that exceptional eigenvalue) will have smaller period and 
we should exclude them from initial data. 



First let us construct the K-matrices which have the eigenvalues in quadratic 
extension GF[yJp\. If h is the primitive element of the GF[y/p], that is 



h = h l + h 2 ^g, h-h* = g, h + h* = 2h h (14) 

then the matrix which has the eigenvalues in GF[^fp\ can be constructed in the form 
of (4) 



A, 



0. 
0. 



1- 
0. 



V -1, Zhg -g, 
because the characteristic equation is 




(mod p) 



(15) 



(\ + g-)(\-h)(\-h*) = 

A 3 - (2/ix - g~)\ 2 - (2h 1 g~ - g)X + 1 = (modp) (16) 

and has two roots in GF[y/p] and one root in GF\p\. The period of the most 
trajectories is equal to 

r p =p 2 -l. (17) 

and is quadratic in p. At the same time the trajectories with the initial vector 
corresponding to eigenvalue A = — g~ are smaller and one should exclude them from 
initial data. It is also easy to see that if we want to construct two-dimensional 
K-matrices with eigenvalues only in GF[*Jp] we face the problem with determinant 
Det A — h • h* — g ^ 1. This observation explains why two-dimensional K-systems, 
like Arnold cat, can not have periodic trajectories of the length p 2 — 1 on any Galois 
sublattice. 

To construct a K-matrix generator with eigenvalues in high field GF[^/p] we will 
use primitive polynomial of degree d over GF[^fp\. The primitive polynomial has 
the form |(J || 

\ d + PiX dr ~ 1 + lhX i - 2 + ... + Pd = (18) 

with coefficients fii, f3 2 , fld over GF[p\. The roots of this characteristic polynomial 
coincide with different powers of a primitive element h (12) of GF[tfp] 



Ai = h, As 
If p d — 1 is not divisible by p, p 2 , ...,p d 



h p , ... ,\ d = h p - 

then all of them are primitive elements 



of GF[jVp\. This is the reason why this polynomial is called "primitive". There are 
two equivalent representations of h: i) in the form of root of the polynomial (18) 



and ii) in the form of corresponding matrix |21| , |20| , |22 
/ 



A d 



0. 
0. 



-Pd 




1 



-fa 



01 



(mod p). (19) 



As we already explained the problem is that the primitive polynomial (18) and the 
corresponding primitive matrix (19) do not have determinant equal to one, because 
ft ^ 1. But this property is incompatible with K-condition (2). The exceptional 
case is only GF[2] . 

Nevertheless one can solve this problem as follows: the last term (3d which is equal 
to the determinant of the primitive matrix coincides with the primitive element g 
of GF\p] (3d = g, therefore if we multiply the primitive polynomial (18) by A + g~ 
we will get the polynomial 



(A + g-)(X d + /5iA d_1 + (3 2 X d ~ 2 + ... + (3 d ) 



A d+i + (A + g-)\ d + (ft + /3ig-)X d - L + ... + 1 = 0. 
to which corresponds a matrix with unit determinant of the form (4) 



/ 0, 1 , 
0, 0, 1 



A 



d+l 



o \ 





(20) 



(mod p). (21) 



0, ...0,1 

V -i, ,-(& + fag-), -iPx + g-)) 

of dimension d+l. The trajectories generated by this matrix will have the period 



r p = p d - 1 (22) 

and we should exclude "dangerous" trajectories corresponding to eigenvalue A = 
— g~. They have the form = (x\, x\(—g~), x\(— g~) d ) and very short period 
(p-l)/2. 

Fascinating result of Niedereiter [jIH] allows to estimate the uniformity of maxi- 
mally long trajectories in terms of Kolmogorov discrepancy Dn 

Dr v 1 



The result is very important because the convergence of the Monte-Carlo simulations 
essentially depends on Dn 



10 



The example of the primitive polynomial on GF[\/7] with c? = 10 is [20| A 10 + 



A 9 + A 8 + 3 = and (20) has the form A 11 
the matrix is 



A 



10 



A 9 - 2A 8 -4A + 1 = therefore 



/ 0, 
0. 



.4 



li 



1 






1 



o \ 





V 



0,0,0, ...,0,0,0,1 
-1,4, 0,..., 0,2, 1,1 J 



(mod 7) 



(23) 



and the trajectories have the period 7 10 — 1 except of two trajectories with the initial 
vectors of the form X<® = (1, 2, 4, 1, 2, 4, 1, 2, 4, 1, 2) and X<® = (3,6,5,3,6,5,3,6,5,3,6). 
The same matrix will have different properties on Galois field GF[p'] where p' ^ p. 



The determination of the set of primes for which a given matrix has the maximal 
period is an unsolved problem |33] . 

Tables of primitive polynomials with large values of d are available 120 
contains tables for d < 101 



particular 37 



in 



for d < 169 and in 39 



In 

for 



d < 1001 with the corresponding periods of order 2 



1000 



4.1 Conclusion 

In this article we advocate two approaches to generate pseoudorandom numbers of 
high quality: i) the first one is based on K-system generators with their exponentially 
unstable trajectories uniformly filling the phase space and ii) on primitive matrices 
acting on a given Galois sublattice with their maximally long trajectories. We 
demonstrate that one can combine these properties in a unique K-matrix which has 
primitive matrix as a submatrix. This construction guarantees that the trajectories 
belong to a K-system and at the same time have maximally large periods. 
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